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Abstract. The authors investigate the solution of a nonlinear reaction- 
diffusion equation connected with nonlinear waves. The equation discussed 
_j , is more general than the one discussed recently by Manne, Hurd, and Kenkre 

>• I (2000). The results are presented in a compact and elegant form in terms 

of Mittag-Leffier functions and generalized Mittag-Leffler functions, which 
Tj- ■ are suitable for numerical computation. The importance of the derived re- 

^ . suits lies in the fact that numerous results on fractional reaction, fractional 

\Q , diffusion, anomalous diffusion problems, and fractional telegraph equations 

O I scattered in the literature can be derived, as special cases, of the results 

investigated in this article. 






c^ 



1 Introduction 



r> I Reaction-diffusion models have found numerous applications in pattern for- 

d ' mation in biology, chemistry, and physics, see SmoUer (1983), Grindrod 

(1991), Gilding and Kersner (2004), and Wilhelmsson and Lazzaro (2001). 
These systems show that diffusion can produce the spontaneous formation 
spatio-temporal patterns. For details, refer to the work of Nicolis and Pri- 
gogine (1977), and Haken (2004). A general model for reactiondiffusion sys- 
tems is discussed by Henry and Wearne (2000, 2002), and Henry, Langlands, 
and Wearne (2005). A piecewise linear approach in connection with the diffu- 
sive processes has been developed by Strier, Zanette, and Wio (1995) which 
leads to analytic results in reaction-diffusion systems. A similar approach 
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was recently used by Manne, Hurd, and Kenkre (2000) to investigate effects 
on the propagation of nonlinear wave fronts. 

The simplest reaction-diffusion models can be described by an equation 

ON d'^N 

-^0—^,F(N), (1) 

where D is the diffusion coefficient and F{N) is a nonlinear function rep- 
resenting reaction kinetics. It is interesting to observe that for F{N) = 
7iV(l — A^), eq.(l) reduces to Fisher-Kolmogorov equation and if we set 
F{N) = '~fN{l — N'^), it reduces to the real Ginsburg-Landau equation. 

A generalization of (1) has been considered by Manne, Hurd, and Kenkre 
(2000) in the form 

d^N ON ^d^N .^^, , 

where ^ indicates the strength of the nonlinearity of the system. In this arti- 
cle, we present a straightforward method for the systematic derivation of the 
solution of nonlinear reaction-diffusion equations connected with nonlinear 
waves, which is more general than the equation (2). The results are derived 
in a closed-form, by the application of Laplace and Fourier transforms, which 
are suitable for numerical computation. The present study is in continuation 
of our investigations reported earlier in the articles (Saxena, Mathai, and 
Haubold, 2002, 2004, 2004a, 2004b, 2005). 



2 Mathematical Prerequisites 

A generalization of the Mittag-Leffler function (Mittag-Leffler 1903, 1905) 

oo n 

^"W = Ef7 —TV «eC,i?e(a)>0, (3) 

^0 r(na + 1) 

was introduced by Wiman (1905) in the general form 

oo ~n 

EaA^) = J2 FT — ^' a,/3eC, Reia) > 0, i?e(/5) > 0. (4) 

^0 r(na + (3) 



The main results of these functions are available in the handbook of 
Erdelyi, Magnus, Oberhettinger, and Tricomi (1955, Section 18.1) and the 
monographs by Dzherbashyan (1966, 1993). Prabhakar (1971) introduced a 
generalization of (4) in the form 

KM = E r( i%V ^1 "'^'^ ^ C;Reia),Re{P),Re{^) > 0, (5) 
n=o y^'^ ' P)\^)- 

where (7)^ is Pochhammer's symbol, defined by 

(7)0 = 1, (7)r = 7(7 + 1)(7 + 2) . . . (7 + r - 1), r = 1, 2, . . . , 7 ^ 0. (6) 

It is an entire function with p = [Re{a)]^^ (Prabhakar, 1971). The so- 
lution of generalized Volterra-type differ-integral equations associated with 
this function as a kernel is derived by Kilbas, Saigo and Saxena (2004). A 
general theory of generalized fractional calculus based on this function has 
been developed by Kilbas, Saigo and Saxena (2004), generalizing the results 
for Riemann-Liouville fractional integrals and derivatives, which form the 
backbone of fractional differ-integral equations. 

For 7 = 1, this function coincides with (4), while for /5 = 7 = 1 with (3) 

Elp = E^,p{z),El,{z) = E4z). (7) 

We also have 

^a,f3;z)= ^F,ia,(3;z) = Tif3)Elf,iz), (8) 

where ^{a,(3;z) is Kummer's confluent hypergeometric function defined in 
Erdelyi, Magnus, Oberhettinger, and Tricomi (1953, p. 248, eq.(l)). Prab- 
hakar (1971, p. 8, eq.(2.5)) has shown that 

/■oo 

/ t^-'e-''El^{ujt'')dt = s-"'{l-us-^)-\ (9) 

where Re{P) > 0,Re{'j) > 0,i?e(s) > and s > |c^|«^(«. 

The Riemann-Liouville fractional integral of order z/ is defined by Miller 
and Ross (1993, p.45) 

oDrfit) = :p^ f\t - ur~^f{u)du, Re{u) > 0. (10) 
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Following Samko, Kilbas, and Marichev (1990, p. 37), we define the frac- 
tional derivative for Re{a) > in the form 

* •^ ^ ^ T{n-a) dr Jo {t - m)°-"+i ' l j > \ > 

where [a] means the integral part of the number a. In particular, if < a < 

„«?/(«) = 1^-^/'^. (12) 

at r(l — a) JO [t — u)°' 

and if q; = n G A^ = {1, 2, . . .} , then 

oD^f{t) = D-fit),iD = d/dt), (13) 

which is the standard derivative of order n. From Erdelyi, Magnus, Ober- 
hettinger, and Tricomi (1954a, 1954b, p. 182), we have 

L{oD;''f{ty,s} = s-''f{s), (14) 

where /(s) is the Laplace transform of f{t), defined by 

L {fit); s} = fis) = / expi-st)fit)dt, Re{s) > 0, (15) 

Jo 

which may be written symbolically, as follows 

f{s) = L {fit); s} or f{t) = L'' {/(s); t} , (16) 

provided that the function f(t) is continuous for t > and of exponential 
order as t ^ oo. The Laplace transform of the fractional derivative is given 
by Oldham and Spanier (1974, p. 134, Eq. (8.1.3)) 

n 

L UD^fit); s} = s"/(s) + Y: s'-' oDr^/m^o. (17) 
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In certain boundary-value problems, the following fractional derivative of 
order a > of a causal function f{t) (that is, f(t) = 0, for t < ), is 
introduced by Caputo (1969) in the form 

1 /■* f'^'^Hr)dT 
D^f(t) = — -/ -^ ^-^—T,m-l<a<m,Re(a)>0,meN,. 

d"^f 
= — — , if a = m, (18) 



where -^^f is the m^^ derivative of /. 

Caputo (1969) has given the Laplace transform of the fractional derivative 

as 

m—l 

L {Dff{t); s} = sTis) - ^ s°-^-V('')(0+), m - 1< a < m, (19) 

where F{s) is the Laplace transform of f{t). 

The above formula is useful in deriving the solution of differ-integral equa- 
tions of fractional order governing certain physical problems of reaction and 
diffusion. We also need the Weyl fractional operator defined by 

where n = [fi] is an integral part of /i > 0. Its Fourier transform is given by 
Metzler and Klafter (2000, p.59, A.ll, 2004) 

F{_^D^JXx)} = itkrnk), (21) 

where we define the Fourier transform by the integral equation 

/oo 
h{x)exp{iqx)dx. (22) 

-co 

Following the convention initiated by Compte (1996), we suppress the imag- 
inary unit in Fourier space by adopting a slightly modified form of the result 
(21) in our investigations (Metzler and Klafter 2000, p. 59, A. 12, 2004) 

F{^^D^JXx)} = -\krnk). (23) 

Now we will establish the following results, which provide the inverse Laplace 
transforms of certain algebraic functions and are directly applicable in the 
analysis of reaction-diffusion systems that follows. 
It will be shown here that 

-^^-j^; <| = i;(-a)'«'"-''E:S_,,„(-«"). (24) 



where Re{a) > 0,Re{f3) > 0,Re{s) > 0, |^^| < 1, and -E^^^(x) is the 
generahzed Mittag-Leffler function defined in eq. (5). 
Proof. We have 



„a-l „a-/3-l 



oo f ^\r 

_ ea-/3-l 



Taking the inverse Laplace transform of the result (25) and using eq. (9), we 
obtain the result (24). The term by term inversion is justified by virtue of 
a theorem given by Doetsch (1956, 22). In a similar manner, the following 
results can be established 

^ ' [s'^ + as'^ + b } 

oo 

= t--'T.i-^yt^''~'^'K:ii-,)ir^i)^ii-bn, (26) 

r=0 

where Re{a) > 0,i?e(/3) > 0,Re{s) > 0, |-^| < 1, and a> (3. 






e-' T.i-^yt^'"'^''EZaMc.-,)ri-bn, (27) 



r=0 
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where Re{a) > 0,Re{/3) > 0,Re{s) > 0,and \-^f^\ < 1. 
We now show that 



(^) ^'' \s2- + as- + b'^\ 



, ^ [(A + a)E,(At") - (/i + a)E„(/it")], (28) 

v/(«' - 4&) 



where a^ — 4b > and Ea{x) is the Mittag-Leffler function defined in eq. 
(3), Rei^a) > 0, Re{s) > 0, and A and /i are the real and distinct roots of the 
quadratic equation, x'^ + ax + b = 0, namely A = |(— a + \/{a^ — 46)) and 



l(-a-^ia^-Ab)). 



Proof. We have 



S2"-1 + as'^-^ 1 ,5^"-^ + as""-^ S^a-l + Qga-1_ 



s2" + as" + 6 A - ;u ^ s° - A s" - /i 

1 (A + a)s"-i (/i + a)s"-i 
~ A-;U s" - A s"^ - jj, ^' 

Taking the inverse Laplace transform of (29) gives 

1j 



(29) 



52" + as"' 



^ [(A + a)E,(At°)-(/i + a)E„(/it")]. 



\ — fi 



This completes the proof of eq. (28). 

In a similar manner, it can be shown that 

(^) ^'' \ 2aja^j} = ^[EaA>^n " K,.(/xt")], (30) 

I s^" + as" + J A — /i 
where Re{a) > 0, Re{s) > 0, and A and n are given along with (28). 



3 Solution of Fractional Reaction-Diffusion 
Equation 

In this section, it is proposed to derive the solution of the fractional- diffusion 
system connected with nonlinear waves governed by the eq. (31). This system 
is a generalized form of the react ion- diffusion equation recently studied by 
Manne, Hurd, and Kenkre (2000). The result is given in the form of the 
following theorem. 
Theorem. Consider the fractional reaction-diffusion equation 

oD^N{x, t)+a oD^N{x, t) = u^_^D2N{x, t) + ^^N{x, t) + ip{x, t), (31) 



X e?R:,t>0,0<a<l,0<(3< 1, 
with initial conditions 

N{x,0) = f{x),ioT{xe^), (32) 

where z/^ is a diffusion coefficient, y? is a constant which describes the nonhn- 
earity in the system, and <f{x,t) is a nonhnear function for reaction kinetics, 
then there holds the following formula for the solution of (31) 

N{x,t) = 5^^-^/ t("-^)" r(A;)ea;j9(-A;x) X 

r=0 -'-oo 

X [EZil-,)r+i i-bn + t"-'^i5;;;L,)(.+,Hi (-for )]rfA; 

^ y-^J^ ^a+(a-f3)r-l ^{k,t - i)exp{-ikx) X 

r=0 ^TT "'0 ^-oo 

X i^r.+(.-,).(-&r)rfA:rfe, (33) 

where a > /9 and E^^^{-) is the generalized Mittag-Leffler function, defined 
in (5) andb = u^\k\^'-^'^. 

Proof. Applying the Laplace transform with respect to the time variable t 
and using the boundary conditions , we find that 

s"N{x, s) - s°"V(a;) + as^N{x, s) - as^-^f{x) 
= u^ -^D'^N{x,s)+eNix,s) + f{x,s). (34) 

If we apply the Fourier transform with respect to the space variable x, it 
yields 

s"iV*(A;, s) - s'^^Tik) + as^N*{k, s) - as^'^f*ik) 
= -u^\k\m*{k,s) + ^'^N*{k,s) + f*{k,s). (35) 

Solving for N*{k, s), it gives 

where b = v'^\k\"' — ^'^. To invert eq. (36), it is convenient to first invert the 
Laplace transform and then the Fourier transform. Inverting the Laplace 



transform with the help of the results (28) and (30), yields 



r=0 

X K%-p)rU-^ndi. (37) 

Finally, the inverse Fourier transform gives the desired solution in the 
form 

°^ (—nY r°° 



r=0 



X K{a~p)(r+i)+ii-bt'")]expi-ikx)dk 

°° (—nV ft , , foo 

+ E V^ / C^^'^-^^'-' / expi-tkx)^*ik, t - X 

„_n 27r JO J-oo 



X i^ri(.-m.(-&r)c^A;rfe (38) 



r=0 

pr- 
''a,a+{a~/3)r 

This completes the proof of the theorem. 



4 Special cases 

When f{x) = S{x), where 6{x) is the Dirac delta function, the theorem 

reduces to the following 

Corollary 1. Consider the fractional reaction-diffusion system 

oD^N{x, t)+a oD^N{x, t) = v^ ^^DlN{x, t) + i''N{x, t) + v?(x, t), (39) 

subject to the initial conditions 

iV(x, 0) = 5{x) for < a < 1 and < /3 < 1, (40) 

where 5{x) is the Dirac delta function. Here ^ is a constant that describes the 
nonlinearity in the system, and ^p{x,t) is a nonlinear function which belongs 



to the reaction kinetics. Then there exists the following eq. for the solution 
of (39), subject to the initial conditions (40) 

r=0 "'-oo '^ 

X Ea{a-l3){r+l)+li—bt°')]dk 

°° ( — nV ft ,, . poo 

+ E ^^ / r^""~'^^~7 v{k,t-o 

r=0 ^^ •^O -^^"^ 

X exp{-tkx)E:;^^l_^^^^^i-bC)dkd^, (41) 

where b = z/^|/c|'' — ^^. 

Now if we set f{x) = S{x),'y = 2,q; is replaced by 2a and P by a, the 
following result is obtained. 
Corollary 2. Consider the following reaction-diffusion system 

d^"N(x,t) d^N(x,t) 2^2^(x,t) ,2.r. N . N 

with the initial conditions 

N{x, 0) = S{x), Nt{x, 0) = 0, < a < 1, (42) 

(p{x,t) is a nonlinear function belonging to the reaction kinetics. Then for 
the solution of (39) subject to the initial conditions (40), there holds the 
formula 



N{x,t) -- 



/oo 
exp{—ikx)x (43) 

-oo 



27ry^(a2 _ 46) 
X {{\ + a)E^{\r)-{fi + a)Eo,{fit'')}dk 

+ Trh'^'' exp{-ikx)v*{k,t-()x 

Zn Jo J-oo 

X [E^Af^C) - E^AfxC)]dkd^] , 
where A and /x are the real and distinct roots of the quadratic equation 

y^ + ay + b = 0, (44) 

given by 

\ = -(-a+ ^{a^ - 4b)] and /x = - (-a - ^ {a^ - 4b)\ , (45) 
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where 6^ = z/^/c^ — ,^^. 

Proof. In order to prove (43), we replace a by 2a and /9 by a, then eq. (36) 

becomes 

7V.ft.) = ^"''t''^°''+^;''-^'. (46) 

Taking the inverse Laplace transform and using the results (26) and (30), 
yields 



N*{k,t) = [{\ + a)E^(\r)-(fi + a)E^(fie)] (47) 

A — /i 

Jo 

where A and /i are given in (45). The application of the inverse Fourier 
transform to the above equation gives the desired result (43). 

Next, if we set ip{x, t) = 0, 7 = 2, replace a by 2a, and P hj a in 
(31), we then obtain the following result , which includes many known results 
on the fractional telegraph equations including the one recently given by 
Orsingher and Beghin (2004). 
Corollary 3. Consider the following reaction-diffusion system 

d^''N{x,t) d''N{x,t) ^d^N{x,t) .^^, , 

with the initial conditions 

N{x, 0) = 5{x), Ntix, 0) = 0, < a < 1. (49) 

Then for the solution of (48), subject to the initial conditions (49), there 
holds the formula 

N(x,t) = , ^ X (50) 

27r^(a2 - 46) 

/- +00 
exp{-ikx) {(A + a)Ea{\t°') - (/i + a)Ea{fit°')} dk 
-00 

where A and ^ are defined in (45), b = v'^k'^ — ^"^ and Ea{x) is the Mittag- 
Leffler function defined by (3). If we set ^^ = 0, then corollary 3 reduces to 
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the result, which states that the reaction-diffusion system 
a2°iV(x,t) d^N{x,t) _ 2<92iV(x,t) 



with the initial conditions 



dt' 



V 



dx^ 



has the solution, given by 

1 



0, < a < 1, 



iV(x,t) 



271^ [a 

oo 



46) 



(51) 
(52) 
(53) 



/oo 
exp{-ikx) {(A + a)E„(At°) - (/i + a)E„(/it")} dfc 
-oo 



where A and fi are defined in (45), b = z/^/c^ and Ea{x) is the Mittag-Leffler 
function defined by (3). Eq. (53) can be rewritten in the form 



N{x,t) 



1 

47r 



exp{—ikx) 



E^{\n + (1 - 




(54) 



where A and /i are defined in (45) and Ea{x) is the Mittag-Leffler function, de- 
fined by (3). The equation (54) represents the solution of the time-fractional 
telegraph equation (51), subject to the initial conditions (52), recently solved 
by Orsingher and Beghin (2004). It may be remarked here that the solution 
as given by Orsingher and Beghin (2004) is in terms of the Fourier transform 
of the solution in the form given below. The Fourier transform of the solution 
of the equations (51) and (52) can be expressed in the form 



N*{x,t) = -{{! + 



4z/2p) 



)E^{\n + (1 



4z/2A;2)) 



)^a(/it^ 



(55) 

where A and ^ are defined in (45) and Ea{x) is the Mittag-Leffier function 
defined by (3). Finally, it is interesting to observe that the solution of various 
fractional-reaction and fractional diffusion and fractional telegraph equations 
scattered in the literature can be derived as special cases of the theorem 
established in this article. 



12 



5 Conclusions 

There is a host of reaction-diffusion equations such as eqs. (1) and (2) that 
allow the formation of wave fronts which maintain their shape despite the 
diffusive element of evolution contrary to linear expectation. Such generation 
of waves plays a particularly important role in spatio-temporal processes 
in physical systems, including astrophysical fusion plasmas (Kulsrud, 2005; 
Wilhelmsson and Lazarro, 2001; Gilding and Kersner 2004). Much work has 
been done for the numerical treatment of such equations. 

The motivation for the research reported in this paper is the derivation of 
analytic closed-form solutions of fractional reaction-diffusion equations (31), 
(39), (42), and (48) that give rise to nonlinear waves in a respective physi- 
cal medium. For this purpose the paper summarizes specific techniques for 
Laplace, Fourier, and Mellin transforms, results for Mittag-Leffler functions, 
as well as the applications of Riemann-Liouville, Weyl, and Caputo fractional 
calculus for tackling fractional reaction-diffusion equations. Closed-form so- 
lutions of the equations are given in terms of Fox's function and their behavior 
for small and large values of the respective parameter is derived. 
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